RIMS 



Correlation Functions of the XXZ model for A < 1 



MichioJimbo a , Kei Miki 6 , Tetsuji Miwa c and Atsushi Nakayashiki d 

a Department of Mathematics, Faculty of Science, 
Kyoto University, Kyoto 606, Japan 

b Department of Mathematical Science, Faculty of Engineering Science, 
Osaka University, Toyonaka, Osaka 560, Japan 

c Research Institute for Mathematical Sciences, 
Kyoto University, Kyoto 606, Japan 

d The Graduate School of Science and Technology, 
Kobe University, Rokkodai, Kobe 657, Japan 



Dedicated to Professor Chen Ning Yang on the occasion of his seventieth birthday 



Abstract. A new approach to the correlation functions is presented for 
the XXZ model in the anti-ferroelectric regime. The method is based on the 
recent realization of the quantum affine symmetry using vertex operators. With 
the aid of a boson representation for the latter, an integral formula is found for 
correlation functions of arbitrary local operators. As a special case it reproduces 
the spontaneous staggered polarization obtained earlier by Baxter. 
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§0. Introduction 

In this article we consider the one-dimensional infinite spin-chain 



Hxxz = - £ (a% +1 a% + ^ +1 tr» + Aa z k+1 a z k ) (0.1) 

k— — oo 

known classically as the XXZ model. We shall limit ourselves strictly to the anti- 
ferroelectric regime A < — 1. Our aim is to find an exact expression for the spin corre- 
lation functions (a? 1 ■ ■ -c^*}, where a] = crf 7 af 7 cr^ and (•} = (vac • |vac) denotes the 
ground state average. Here the ground state |vac) means one of the two ground states 
in the anti- ferroelectric regime. In Baxter's paper [1] they are denoted by |±). The 
present paper is based entirely on the framework of the recent work [2] , which discribes 
the realization of the quantum affine symmetry for (0.1) using the g-deformed vertex 
operators. Let us recall below the contents of [2] which are relevant to the subsequent 
discussions. 

The Hamiltonian (0.1) is an operator acting on the infinite tensor product V® 00 = 
■ ■ ■ V®V®V® • • • of the two dimensional space V = C 2 . The space V® 00 admits also an 
action of the quantum affine algebra U — U q (sl2) via the iterated coproduct. A naive 
computation shows that the algebra U provides an exact symmetry for (0.1). Namely 
if A = (q + q~ 1 )/2, then [Hxxz,U'] = where U' = U' v (sl 2 ) denotes the subalgebra 
'without the grading operator d 1 [2], while d plays the role of the boost operator. However 
the actions of (0.1) and U are both defined only formally, and the issue is how to extract 
the theory free from the difficulties of divergence. 

The basic idea in [2] is to replace the formal object V® 00 by the level [/-module 

Tx.,» = V(\)®V(fi)* a ~ Hom(V(/i), V(X)), (0.2) 

where V(A) (A = A , Ai) denotes the level 1 highest weight [/-module, and V(X)* a signi- 
fies its dual (the superfix indicates that the [/-module structure is given via the antipode 
a). The choice A = /i (resp. A 7^ fx) is responsible for the even (resp. odd) particle 
sector. V(Aq) (resp. V'(Ai)) means that we are working in the boundary condition 
CT 2fe = = - 1 ( res P- a 2k = ~ l ^2k+i = !) for k » °> and ^( A o)* Q (resp. 

V(A 1 )* a ) means the same thing for k « 0. The tensor product should be completed 
in the g-adic sense to allow for infinite sums (see [2] for a precise treatment). To make 
contact with the naive picture of V® 00 , one utilizes the embedding of V(\) into the half 
infinite tensor product • • • V®V®V [3] . This is supplied by iterating the vertex operators 

$f : V(A) — >V{p) ® V. (0.3) 

(It is conjectured [2] that there is a unique normalization of (0.3) which makes the infinite 
iteration convergent.) Similarly V{fi)* a embeds to the other half infinite tensor product 
V®V®V® ■ ■ •, giving altogether the embedding 

— ► ■■■V®V®V® (0.4) 

Eq. (0.4) provides the principle of interpreting the notions defined in the picture V®°° 
by pulling them back to T\^. The translation operator T is the first such example. (The 
Hamiltonian (0.1) itself is defined in terms of T and the grading operator d; see [2].) In 
this paper we follow the same principle to formulate the local operators as acting on 
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In the language of (0.2) the ground state (vacuum) vector is given by the identity 
element of = Hom(V(A), ^(A)). The inner product of two vectors f,g E is 
given as the trace 

(f\o) = try(A)(g ; 2 7f ) (o.5) 



where p = Aq + Ai (we have changed the normalization from [2], see Sect. 1 below). Thus 
the correlation functions (vacla^ 1 ■ ■ ■ oj^ |vac) can be expressed as the trace of products 
of the vertex operators (0.3). 

To perform the evaluation of these traces, we invoke the bosonization method. The 
realization of the q-deformed currents on level one modules was done in [4] using (ordi- 
nary) bosons. In the same spirit we derive the formulas for the vertex operators (0.3) 
in terms of bosons. This leads to an explicit formula for the correlators in terms of 
certain integrals of meromorphic functions. We verify that in the simplest case of the 
one-point function (cr|) this formula reproduces the known result for the spontaneous 
staggered polarization due to Baxter [1]. Jacques Perk noted a strong similarity between 
our formula and the formula for the Ising model spin-spin correlation functions given in 
[5,6]. 

The bosonization of the vertex operators also enables us to compute the n-point 
functions discussed in Sect. 6.8 of [2]. They are the matrix elements (not the trace) 
of the product of the vertex operators with respect to the highest weight vectors. The 
conjectural formula (6.39) of [2] is thus proved. In this paper, we do not go into details 
on this matter. 

The plan of the paper is as follows. In Sect. 1 we formulate the local operators 
and their correlators using the scheme mentioned above. In our algebraic formulation 
the algebra U, intertwiners, etc. are a priori defined over the base field F = Q(q r ) (q 
an indeterminate, see the remark at the end of Sect. 2); we shall see however that the 
resulting formulas are meaningful for complex values of q with \q\ < 1. In Sect. 2 we 
outline the bosonization of the vertex operators. The basic ingredients are the Drinfeld 
realization of the algebra U' and the result of [4], which we recall briefly. Unlike the 
classical case (q = 1) only one of the components <&-(z) of (0.3) has the exponential 
form, while the other one does not but is given as the g-commutator of $-(z) with a 
generator fx of U. Sect. 3 is devoted to the formula for the correlators. 
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§1. Vacuum expectation values of local operators 

Let us recall some notations of [2] . We denote by V — Fv + © Fu_ the 2-dimensional vec- 
tor space on which the Pauli matrices a x 7 a y , a z act. We consider V as the 2-dimensional 

U' q (sl(2)) -module. The AXZ-Hamiltonian formally acts on the infinite tensor product 

V®°° of V. We label the components of the tensor product by integers k € Z from right 
(k — > — oo) to left (+oo <— k). In our mathematical scheme, we replace the semi-infinite 

tensor product of the components k > 1 by the irreducible highest weight C/ 9 (s[(2))- 
modulc V(Aj) (i = 0, 1) Let us define the action of the local operators on V(Ai). Let 
L e End( y <g) ■ ■ ■ © y ). The operator L naturally acts on the tensor product of the com- 

n 

ponents n > k > 1 of We want to interprete this action as one in End(V(Aj)). 

For this purpose, we use the vertex operators [2] : they are the intertwiners 

$f(z) : V(X) -» V(ji) ® $f(z){v) = ® v+ + &£_{z){v) ® w_ (1.1) 

$^(z) : K(A) ® V z -> y(/i), *£ v (s)(« ® v±) = &Z v ±(z)(v), 

where A = A a ,p = Ai or A = Ai,^ = A , and K = V ® F[z, z' 1 ] signifies the U q (sl(2))- 
module associated to V with the spectral parameter z (see [2], Sect. 6). Precisely speak- 
ing, we consider the vertex operators as the generating functions of their Fourier com- 
ponents in terms of the spectral parameter z. They are normalized as 

*fc V (z){\u Ao )) = \u Al )®V- + ---, $ A f r (z)(\u Al )) = \u Ao )®v + + ---, (1.2) 

Here (ma*) denotes the highest weight vector of V(Aj). We have used the standard 
notation (z^p)^ = I1°!L (1 - zp>). 

Set A n = A„_2 for simplicity of notation. Given an L as above, wc define the 
operator C = p£,..., Zl (L) G End(y(A,)) by 

(„2. „&\n 

'■^^•"••♦to' 1 - 1 

In this paper, we use the following convention (different from that of [2]) for the 
invariant bilinear form on P = ZAo © ZAi © ZS: 

(A , A ) = 0, (A , ai) = 0, (A , 5) = 1, 
(ai,ai) = 2,(ai,(5) = 0,(5,5) = 0. 

Note that Ai = Ao + oji/2, 5 = ao + ot\. We set p = Ao + Ai as usual, and also a — ai 
for simplicity. We identify P* = Zh © Z/ii © Zc? as a subset of P via ( , ). We have 
p = 2rf + ai/2. 
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The vacuum expectation value (L) z l },..., Zl of the local operator L is given by 

_ trv( Ai) ( g - 2 ^„„,, 1 (L)) 
tr y(A,)(9 p ) 

This is a consequence of the formula for the invariant inner product in the space V^A^) ® 
V(Ai)* a [2]. For the XXZ model correlator we specialize the spectral parameters to 
z\ = ■ ■ ■ = z n . We expect that the formula unspecialized would give the equal-row verti- 
cal arrow correlator for the inhomogeneous six-vertex model with z k being the trigono- 
metric spectral parameter of the fc-th vertical line. Here we use the usual language of 
vertex models on the 2-dimcnsional square lattice, in which the fluctuation variables are 
described as arrows sitting on vertical or horizontal edges. Note the following selection 

rule, which is specific to the six vertex model: (L)^,..., Zl — for L = E £ > n£n <g> • • • <g> £7 £ < 

(Eij is a matrix unit) such that e\ + • • • + e n ^ s[ + ■ ■ ■ + s' n . 



£, £1 



§2. Bosonization 

Let us first recall Drinfeld's realization of the quantum affine algebra U 1 = U' q (s[(2)} [7]. 

It is an associative algebra generated by the letters {x^ | k G Z}, {a/ | I G Z^o}, 7 ±1//2 
and K , satisfying the following defining relations. 

7 ±1 / 2 G the center of the algebra, 

[a k ,ai} = S k+lfi -[2k] - - , [a k ,K]=0, 

Kx^K- 1 = q ±2 x^ 

[a k ,x±]=± 1 -[2k] 1 ^ x ± +l , 

x t+l x t ~ 9 ±2x t x k+l = Q ±2x k x l+1 ~ x t+l x k' 

where [n] — (q n — q~ n )/(q — q^ 1 ) and {ip r ,tp s \ r G Z> ,s G Z< } are related to 
{ai | I G Z #0 } by 

OO OO 

^2^ k z- k =Kexp^(q-q- 1 )^2a k z- k ^, 

k=0 k=l 

oo oo 

Y V>-kZ k = R- 1 exp(-(g - q- 1 ) Y a_ fc z fe ) . 

fc=0 k=l 

The standard Chevalley generators {e^,/,,^} are given by the identification 
to = 7-^~\ ti = K, ei = x£, fi = Xq, eoh = x± , if/Vo = xt^ 
We use the coproduct 

A(e») = ei <g> 1 + U <g> e», 

A(/ i ) = / i ®*r 1 + l®/i, 
A(ti) = t< (8> *<- 
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Next recall the construction of the level one irreducible integrablc highest weight 
representations of U' in terms of bosons [4]. Let P = Z^, Q — Za be the weight /root 
lattice of sl(2), and let F[P], F[Q] be their group algebras. The basis elements of F[P] 
are written multiplicatively as e na (n G |Z). As an F[Q]-module, F[P] = F[P] © 
F[P]i where F[P]i = F[Q]e la / 2 . Introduce a [/'-module structure on the space W — 
F[a-i, a_2, • • •] <8> ^[-P] in the following way. First let {a,k}, e 13 and d a act on W as 

afc = the left multiplication by ® 1 for fc < 0, 

= [a/j, •] <g> 1 for fc > 0, 
^(/«e*) = /®e /3l+/32 , 



d a (f®eP) = {a,l3)f®eP. 



We let also 



A = 1 <g) ^ , 7 = g <g> id. 



The actions of the generators {x„} are given through the generating functions X ± (z) — 
E„ e z^ z_ " _1 as follows. 

OO OO 

X+(z) = cxp(]T <T n/ V) exp(- £ jV»/ 2 *-")e<"s a <», 

n— 1 L J n— 1 ^ ^ 

oo oo 

X-(,)= e xp(-^^g«/V)exp(^^ 3 

n=l L ' n=l L J 



/ 2 z-")c- Q z- a - 



With these actions W becomes a [/'-module. The submodules F[a_i, a_2, • • •] <8> -F[-P]i 
are isomorphic to the irreducible highest weight modules F(Aj) with the highest weight 
vectors ma = 1 ® 1 and = 1 ® e Q / 2 . The grading operator d is introduced by 

-d(a- • • • a!X e") = (£ n^- + <M - (a- • • • a- e* 3 ). 

i=i 

We remark that the trace of p _d A + (z)A~(z) on F(A ) is calculated in [8]. We shall 
describe the vertex operators in terms of the representation constructed above, and then 
compute a similar trace for them (see Sect. 3). 

Let the components $±(z) = be defined as in (1.1). From the condition 

that it intertwines the action of Xq we find 

*+(*) = [*-(*),*o]*> (2- 1 ) 

where [X, Y] q = XY — qYX. The intertwining relations with the Chevalley generators 

determine $>±(z) uniquely, but it is not easy to get the expression of the vertex operators 
from them only. We better use the relations with Drinfeld's generators a n . We lack the 
formulas for the coproduct of them, in general. However we have the following partial 
information which will suffice for our purpose [9]. 
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fc-i 



A(x+) = x+ ® 7 fe + 7 2fc ^ ® x+ + ^ 7 ( fc+3i >/Vfc-i ® 7 fc -*x+ mod iV_ ® iV 2 , 



»=o 
i-i 



A(a;+ ) = x± t ® 7"' + ® x+_ x + l (l ~ i)/2 <fi-l+i ® J~ l+i xti mod N_ ® JV 2 , 

i=l 

A(aj) = a; ® 7 i/2 + 7 3 ' /2 ® a ( mod iV_ (g) 
A(a_ ; ) = a_; ® 7 _3 ' /2 + 7~' /2 ® a_; mod iV_ ® iV+. 

Here iV± and N± are left F^y ± ,ij) r , p s \r, — s 6 Z>o]-modules generated by {x^\k G Z} 
and {x^x^|m,n G Z} respectively. 

By using Proposition 2.1 and noting that iV±i>± = 0, N + v_ C F[z, we get 

the exact relations 

[a t ,i_(z)]=g a / 2 |lz fc i_(2) fc>0, 
k 

[a_ fc ,$_(z)]= g - 5fc / 2 ^z- fc d_(z) fc>0, 
[$_(z),X+H] =0, t^-(z)t^=q^-(z). 

These conditions along with the normalization (1.2) determine the form of $-(2) com- 
pletely. Explicitly we have, on V(Ai) (i = 0, 1), 

00 00 

$_(*) = CX P (]T ^, 7 "/ 2 z") exp(- V ^- q -^/2 z -ny/2 { _ q 3 z) (d a+ i)/2^ y 2) 

[znj [2tiJ 

n=l L J n=l L J 

It can be checked directly that the operators given by (2.1,2.2) enjoy the correct inter- 
twining properties. 

Let us make a comment on the base field F — Q(q). In the above expression of 
the bosonization of the vertex operators, the square root of q appears. But, it can 
be absorbed in the boson oscillators if we change the definition of their commutation 
relations. Note also that the combination (d a + i)/2 in the power of — q 3 z produces only 
an integer power acting on V(Aj). In any event, the appearance of q 1 ' 2 is only superficial, 
and our theory is free from the ambiguity in the choice of the square root. 
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§3. Calculation of correlators 

As explained in Sect.l, for the evaluation of the correlators it is necessary to calculate 
the trace of the product of the vertex operators. Using the boson realization described in 
Sect. 2, we shall treat the trace of the form tr(x~ d y a O) where x, y are complex parameters 
with \x\ < 1. Since q~ 2p = q~ 4d ~ a , the choice x = q A , y = q^ 1 will be relevant to the 
correlators. The calculation is simplified by the technique of Clavelli and Shapiro ([10], 
Appendix C). Their prescription is as follows. Introduce a copy of bosons {b n } satisfying 
[a m ,6„] = and the same commutation relations as the a n . Let 

For a linear operator O = 0({a n }) on the Fock space T a = F[a-i,a~2, ■ • •], let O = 
0({a n }) be the operator on T a ® Tb {Fb = F[b-i,b-2, ■ • •]) obtained by substituting a n 
for a n . We have then 



n~=i(i-^)' 

where d — — ^2 n>0 r^Tp^l B -« B n an d (0|O|0) denotes the usual expectation value with 

respect to the Fock vacuum |0) = 1 ® 1 € T a ® Tb, (0|0) = 1. 

In order to apply their method to our case, we need the following. Set 

f(z) = (zx; x) QO (q 2 x/z; x) x , g{z) = JJ j^nl^j^ - 

Let us denote the operators X~(z) and $±(z) with a n substituted for a n by J(z) and 
(f>±(z). Then the explicit expressions of J(z) and <f>-(z) (after the normal ordering) are 

00 q n/2 



J(z) = /(l) exp(- V \ T {z n a. n - z ~ Uh -n)) 

q n l 2 z-"a n - (xz) n b n . _ a _ d 

x exp > -—- e z Ua 

v ^ n l-x n ' 



n=l 



,.. n . q^l 2 z n a- n -q^ n l 2 z- n b^ 

n— 1 L J 

x eM-E" [2 %{l-*») V /a (-g 3 *)< a ^ /2 on V(M). 

They satisfy the following relations, 

^(6) I) (-f )/(|)/(|). 



0-(2i)</>-(2 2 ) =: 0_(zi)0_(z 2 ) 



"("Vf— 



(i- 9 2 w -i)/H' 



(l-tu)/M 



( ^ + (z) = ( (Z - (Z - 1 ) & ^ : j(^_(z) 



g2<| w |<l 27Ti' ■ (l-w)(l-g 2 W r )f{w)' 



Correlation Functions of the XXZ model for A < — 1 



9 



Here w — £/ ' q 2 z and : • • • : denotes the normal ordering with respect to the bosons a n 
and b n (We do not normal order the operators e na and d a ). In the above equations, all 
factors of the form 1/(1 — z) should be understood as E^Lo z ™- This fact specifies the 
contour for the expression Using the formulas above, (1.3) is evaluated as follows. 
Set 



p £ e r:;.f{z n ,---,z 1 \x, y \i) 



(<z 2 ;<? 4 )? 



£i i „ „ i „, „ , i -\ y± ; 1 Joe 

4. „4Vi 



(<z 4 ;<? 4 ) 



oc 



x tr y(A0 (x-y< :+iy£ , (zQ ■ ■ ■ ^;-^(zn)^:;! ie „(zn) ■ ■ ■ g(j 

tr v(Ai) (x-V) 

We have 

for L = E e / nEn ® • • • <g) P e < i£l (see the remark at the end of Sect. 1). 
Introduce the following notations 

A = {oi, • ■ ■ , o a } = {j | e ■ = -1}, S = {6i, ■■■M} = {j I Ej = +1}, 
(s + t = n, at < dj, bi < bj for i < j), 
h(z) = (q 2 z; x) OQ (xz; x) OQ (q 2 2T 1 ; x) 00 {xz~ 1 ; x)^. 

We prepare the integration variables £ a (a S A), (& (6 6 P) and set r]j = £, aj (1 < j < s), 
= C&„+i- 3 (« < i <n),fj = Ylj T}j and z = f]j Zj- Then we have 

^"'.■.?; i (^,---,^l^yh) 



) * <z E ae ^+E 6eB "-™(«+ 1 )/ 2 TT / ^ TT / ^ 

ilt2^ a -z )l. 
z j - 9 2 £a TT TT & ? 2 ^' TT Vk — Vj 



ii ii n n ^z^n„, .,,„. 

a£Aa<j<n J Sa beBb<j<n S J j<fe ' y " 

Uj,k h (Vj/zk) (x;x) 00 tr v(Ai) {x- d y a ) 
Note that the last factor of the above equation can be rewritten into 

'zy{-(yz/f)) 2 x 1+i \x 2 ) oo {-{fj/yz) 2 x 1 - i \x' 2 ) O0 



77/ (-y 2 a; 1 + i ;a; 2 ) 0O (-y 2 x 1 'ji^ 

The contours of integration should be chosen as follows. Both C a and Cb are anti- 
clockwise, and the z% (1 < i < n) lie inside of C a and outside of Cb- Other relevant poles 
with small multipliers q, x (e.g., q 2 Zi) are inside, and those with large multipliers are 
outside of the contours C a and Cb- By specializing (3.2) to the case x = q A and y = q -1 , 
we obtain the integral formula for the correlators. The case n = 1 gives the spontaneous 
staggered polarization P = P+ (1 | 1) -PT(1 1 1) [1]. Note that P/(l | i) = PT e e (1 | 1 — *). 
By using (3.2) for P+(l | 1) and P+(l | 0) we get P = (<? 2 ; g 2 )^/(-g 2 ; g 2 )^. 
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